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tronic  devices  and  nonlinear  optics.  The  enhanced  photoluminescence  effect  in  monolayer 
M0S2  was  discovered  and,  as  a  strong  tool,  was  employed  for  strain  and  defect  analysis  in 
MoS2.  Recently  large-size  monolayer  M0S2  has  been  produced  by  chemical  vapour 
deposition,  but  has  not  yet  been  fully  explored.  Here  we  systematically  characterize  chemical 
vapour  deposit!  on -grown  M0S2  by  photoluminescence  spectroscopy  and  mapping  and 
demonstrate  non-uniform  strain  in  single-crystalline  monolayer  MoS2  and  strain-induced 
bandgap  engineering.  We  also  evaluate  the  effective  strain  transferred  from  polymer 
substrates  to  M0S2  by  three-dimensional  finite  element  analysis.  Furthermore,  our  work 
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High  quality  and  single  layered  molybdenum  disulfide 
(M0S2)  has  been  synthesized  by  chemical  vapour 
deposition  (CVD)  over  the  past  2  years ^  boosting  its 
promise  in  applications  involving  hydrogen  evolution  reaction^®, 
supercapacitors ultra  low  power  consumption  field  effect 
transistors  and  phototransistor^^’^^,  spintronic  devices^^, 
nonlinear  optics^^  and  biosensors^^.  CVD  grown  M0S2  atomic 
layers  are  also  superior  to  samples  prepared  by  mechanical 
exfoliation  and  chemical  isolation  methods,  as  continuous 
monolayer  M0S2  (hereafter  IH  M0S2)  can  currently  be  grown 
at  the  centimetre  length  scale^’^^’^^’^^.  However,  concurrently 
with  the  large  scale  growth  of  monolayer  M0S2,  variations  in 
crystal  geometries  and  atomic  structures  have  been  predicted  and 
observed^^.  Geometrically,  M0S2  triangles,  hexagons,  six  pointed 
stars,  irregular  shapes  and  continuous  films  can  be  synthesized, 
depending  on  the  precise  growth  conditions  such  as  pressure, 
temperature,  substrates,  precursors,  cooling  rate  and  so 
on3  5,18,20  Structurally,  like  graphene,  a  wide  variety  of  grain 
boundaries  (GBs)  has  been  observed  in  polycrystalline  CVD 
grown  M0S2  (refs  3,4,21).  Thus,  it  is  important  to  understand 
how  these  variations  modify  the  physical  properties  of  M0S2. 

It  is  well  known  that  M0S2  has  an  indirect  to  direct  bandgap 
transition  as  it  is  thinned  down  from  multi  layer  to  monolayer, 
producing  a  significant  enhancement  of  photoluminescence  (PL) 
quantum  yield  as  a  result  of  the  direct  bandgap  excitonic 
transition^^’^3  However,  the  intensities  and  peak  positions  of 
such  PL  signals  closely  correlate  with  many  other  parameters 


such  as  crystal  structure,  local  doping  and  strain  within  the 
sample  and  require  more  detailed  examination.  CVD  grown 
M0S2  with  the  richness  of  its  intrinsic  and  extrinsic  defects  show 
significant  heterogeneity  in  PL  characteristics.  Consequently,  a 
systematic  examination  of  PL  characteristics  in  CVD  grown 
M0S2  could  provide  information  to  help  better  understand  the 
growth  mechanisms,  crystal  structure  and  band  structure  of  these 
monolayers. 

In  this  work  we  examine  the  PL  characteristics  of  CVD  grown 
M0S2  monolayers  with  a  focus  on  the  role  of  growth  induced 
non  uniform  strain  and  GBs.  We  examine  the  PL  intensity  and 
peak  position  variations  observed  in  different  samples  and 
explore  the  possible  contributions  to  such  changes.  We  find  that 
non  uniform  tensile  strain  is  intrinsic  to  CVD  grown  M0S2  on 
Si02  and  such  strain  variation  can  be  identified  by  PL  mapping. 
Such  strain  is  induced  through  interactions  with  the  Si02 
substrate  and  can  be  released  after  transferring  M0S2  to  other 
substrates.  On  the  other  hand,  controlled  tensile  strain  could  be 
purposely  introduced  to  monolayer  M0S2  by  stretching  the 
deformable  substrates  it  sits  on,  therefore  providing  an  effective 
method  for  bandgap  engineering.  To  understand  the  role  of 
substrate  properties  on  effective  strain  transfer  and  corresponding 
band  structure  changes  in  M0S2,  three  dimensional  (3D)  finite 
element  analysis  (PEA)  and  DFT  calculations  were  performed. 
These  results  indicate  that  substrates  with  high  Youngs  modulus 
can  dramatically  reduce  the  strain  loss  during  tensile  experiments, 
while  the  commonly  used  polydimethylsiloxane  (PDMS)  can  only 


Figure  1  |  Geometries  and  atomic  structures  of  CVD-grown  monolayer  M0S2.  (a-c)  Optical  images  of  (a)  single-crystal  M0S2  triangle, 

(b)  polycrystalline  l\/loS2  with  nucleation  sites  located  at  the  centre  and  (c)  polycrystalline  six-pointed  l\/loS2  star.  Scale  bar,  5, 10  and  5  pm,  respectively, 
(d)  Optical  image  of  a  M0S2  film  with  second-layer  M0S2  stripes  grown  along  the  GBs,  indicated  by  the  white  arrows.  Scale  bar,  5  pm.  (e)  Annular 
dark-field  scanning  TEM  image  of  defect-free  l\/loS2.  Scale  bar,  0.5  nm.  (f)  Grain  boundary  in  M0S2  with  a  rotation  angle  of  ~17°.  The  grain  boundary 
is  highlighted  by  the  white  dash  line.  Scale  bar,  1  nm. 
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transfer  ~  10%  tensile  strain  to  M0S2.  Finally,  examination  of  line 
defects  in  monolayer  M0S2  samples  shows  that  GBs  can  modulate 
the  PL  response  of  the  material  in  ways  different  from  what  has 
been  previously  proposed.  This  not  only  gives  insight  into 
mechanisms  of  PL  modulation  at  the  GBs  but  also  exemplifies  a 
technique  superior  to  dark  field  imaging  in  transmission  electron 
microscopy  (TEM)  for  their  non  destructive  visualization  in 
M0S2  monolayers. 

Results 

Sources  of  PL  heterogeneity  in  CVD-grown  M0S2.  To  explore 
the  heterogeneity  in  CVD  grown  M0S2  atomic  layers,  we 
synthesize  M0S2  single  crystal  and  polycrystalline  films  using 
M0O3  nanoribbons  and  sublimated  sulfur  as  the  sources  for  Mo 
and  S,  respectively,  and  react  at  a  temperature  of  ~850°C 
followed  by  a  fast  cooling  phase  (see  Methods)^.  The  growth  of 
large  area  films  has  been  proven  to  be  a  process  of  nucleation, 
growth  and  coalescence  of  single  crystalline  triangular  domains. 
More  importantly,  the  formation  of  these  films  is  nucleation 
limited,  and  using  pristine  silicon  dioxide  substrates,  one  can 
truncate  the  process  to  form  isolated  single  crystals.  Our 
experiments  show  that  the  configuration  and  position  of  the 
M0O3  relative  to  the  designated  growth  substrate  is  critical  to  the 
shape  and  structure  of  the  as  grown  crystals.  As  a  result,  one  can 
promote  the  growth  of  polycrystalline  flakes  and  continuous  films 
with  unique  geometries  suitable  for  comprehensive  study  of  PL  in 
M0S2  (Fig.  la  d). 

Some  general  features  of  the  M0S2  flakes  and  films  grown  using 
this  method  are  shown  in  Fig.  1,  such  as  single  crystal  triangular 
M0S2  flake  (Fig.  la),  irregular  M0S2  shapes  (Fig.  lb),  symmetric 
six  pointed  M0S2  star  (Fig.  Ic)  and  polycrystalline  M0S2  films 
(Fig.  Id).  The  bright  areas  located  at  the  centre  of  M0S2  layers 
(Fig.  lb)  are  the  original  nucleation  sites  for  the  crystals.  The 
quality  of  the  CVD  grown  M0S2  is  comparable  to  that  of  the 
exfoliated  samples,  which  has  been  confirmed  by  multiple 
characterization  techniques^.  The  perfect  lattice  of  monolayer 
M0S2  (Fig.  le)  can  be  seen  from  the  annular  dark  field  scanning 
TEM  imaging.  The  brighter  atomic  sites  are  Mo  atoms,  while  the 
weaker  ones  are  double  S  columns.  GBs  are  also  found  in  the 
CVD  grown  M0S2  (Fig.  If)  with  many  variations  as  discussed  in 
our  previous  work^’^L  For  large  area  films,  some  GBs  are 
observed  from  the  optical  images  where  a  second  layer  M0S2 
tends  to  grow  along  the  GBs,  resulting  in  dark  purple  stripes 
indicated  by  the  white  arrows  (Fig.  Id).  Typical  Raman  spectra 
and  intensity  mapping  from  these  mono  /bilayer  topographies 
are  shown  in  the  Supplementary  Fig.  1.  These  features  indicate 
that  GBs  are  preferred  sites  for  nucleation  of  second  layers. 
However,  these  optically  identifiable  features  do  not  fully  account 
for  all  the  GBs  present  in  the  samples  (Fig.  lb,c).  PL 
characterization  of  these  typical  crystals  can  provide  more 
information  about  distinctive  features  of  different  growth 
modes  due  to  the  enhanced  intensity  of  monolayer  M0S2  (ref.  22). 

Native  strain  variations  in  CVD-grown  M0S2  monolayers. 

Figure  2  presents  high  resolution  Raman  and  PL  mapping  of  a 
M0S2  triangular  single  crystal,  acquired  with  long  integration 
times  to  observe  subtle  shifts  in  the  signal.  Raman  mapping  of  the 
E^2g  and  Aig  peaks  shows  a  uniform  intensity  (with  <  5%  var 
iation)  over  the  triangle  (Fig.  2a,b).  However,  significant  varia 
tions  are  seen  in  the  PL  intensity  mapping  (Fig.  2c)  and  it  can  be 
seen  that  the  edges  and  some  regions  within  the  triangles  have  up 
to  a  ~200%  higher  intensity.  Raman  and  PL  peak  position 
mapping  provide  images  with  better  contrast.  Clear  lines  related 
to  the  shift  in  peak  positions  ( ~  1.0  cm  “  ^  hardening  in  E^2g  and 
0.4  cm  “  ^  in  Aig)  can  be  observed  (Fig.  2d,e).  More  interestingly. 


the  best  contrast  is  seen  in  the  PL  peak  position  mapping  (blue 
lines  in  Fig.  2f)  along  and  perpendicular  to  the  edges,  corre 
sponding  to  a  blue  shift  of  ~8meV  (Fig.  2g).  We  propose  that 
these  characteristics  may  arise  from  the  growth  process  that  ends 
with  a  fast  cooling  phase.  Since  the  coefficient  of  thermal 
expansion  of  silica  is  roughly  six  times  smaller  than  M0S2 
(refs  24,25),  a  fast  cooling  process  from  the  growth  temperature 
(~850°C)  could  create  a  significant  contraction  mismatch 
(Supplementary  Discussion).  As  a  result,  the  monolayer  M0S2 
triangles  will  be  in  tensile  strain.  Local  strain  variation  regions 
extended  from  the  edges  vertically  are  towards  the  centre 
(Fig.  2d  f).  In  other  words,  after  cooling,  the  CVD  M0S2 
monolayer  triangles  are  under  a  tensile  strain  but  contains 
localized  regions  where  this  strain  has  been  partially  released  (less 
strained)  possibly  due  to  combined  effects  of  geometry  and 
inevitable  variations  in  interactions  with  the  substrate.  To  esti 
mate  this  CVD  growth  induced  strain  variations,  density  func 
tional  theory  (DFT)  calculations  are  performed  and  a  relationship 
between  the  PL  peak  position  shift  and  the  strain  variations  in  the 
sample  is  established  (Fig.  2h).  The  corresponding  difference  in 
strain  is  estimated  to  be  ~0.2%.  Although  the  DFT  simulations 
do  not  account  for  the  excitonic  effects,  the  relative  change  of  the 
calculated  bandgap  with  respect  to  uniaxial  strain  (  —  44  meV  per 
%  strain)  matches  quite  well  with  previous  experiments  on  the  PL 
peak  position  shift  with  strain  (  —  45  meV  per  %  strain^^  and 
—  36  meV  per  %  strain^^.  Furthermore,  as  we  mentioned,  due  to 
the  fast  cooling,  a  global  tensile  strain  is  found  in  M0S2.  To 
support  our  conjecture,  we  transfer  M0S2  to  a  new  Si/Si02 
substrate  and  compare  the  PL  mapping.  As  expected,  the  tensile 
strain  in  the  film  as  well  as  the  local  variations  are  released  after 
the  transfer  process  (Supplementary  Figs  2  and  3).  The  global 
tensile  strain  is  estimated  to  be  ~  1%.  This  information  provides 
insight  into  the  role  of  thermal  property  mismatch  between  the 
substrate  and  M0S2  in  the  CVD  process.  Our  observations  suggest 
that  PL  studies  can  be  effectively  used  as  a  strain  mapping  tool 
and  a  tool  to  evaluate  strain  induced  bandgap  engineering  if  an 
independent  and  robust  link  between  the  two  parameters  can  be 
established.  To  further  explore  the  role  of  substrate  properties  in 
M0S2  strain  engineering,  a  well  controlled  experiment  is 
performed.  This  experiment  will  also  correlate  strain  with 
optical  and  band  structure  properties  of  monolayer  M0S2. 

Bandgap  engineering  and  effective  strain  transfer  for  M0S2. 

Substrate  induced  straining  of  2D  materials  for  band  structure 
engineering  has  been  effectively  explored  in  previous  studies^^’^^. 
However,  these  experiments  have  only  been  performed  on 
exfoliated  samples  and  a  comprehensive  study  of  the  role  of 
substrate  and  homogeneity  on  the  transferred  strain  has  not  been 
explored.  Most  recently,  the  strain  induced  bandgap  changes  in 
monolayer  M0S2,  using  polymer  substrates,  has  been 
demonstrated^^’^^  However,  because  of  the  weak  adhesion 
between  M0S2  and  substrates,  strain  cannot  be  completely 
transferred  to  M0S2  from  the  substrate,  as  exemplified  in 
graphene  strain  engineering  studies^  Understanding  how 
much  strain  is  effectively  transferred  to  the  M0S2  samples 
becomes  critical  to  bridge  experimental  data  and  theoretical 
calculations.  In  our  work,  monolayer  M0S2  triangles  have  been 
transferred  from  Si02  to  PDMS  (~  0.5  mm  in  thickness;  Fig.  3a). 
As  we  discussed  in  Fig.  1  and  Supplementary  Fig.  2,  through  the 
transfer,  the  thermal  mismatch  induced  tensile  strain  in  M0S2 
and  related  local  strain  variations  are  released.  Uniaxial  force  is 
then  applied  to  the  sample  in  the  direction  indicated  by  the 
arrows  in  Fig.  3a,c  and  was  directly  observed  since  the  M0S2  on 
PDMS  is  visible  under  optical  microscopy  (Fig.  3b).  Deformation 
of  M0S2  flakes  is  clearly  seen  after  the  application  of  force  on  the 
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Figure  2  |  Intriiisk  tensile  strain  In  single-crystal  CVD  MoS^.  (a-c)  Raman  and  PL  intensity  mapping  at  the  Raman  vibrating  modes  and  the 

strongest  PL  peak  from  the  Ae)(citorv  respectively.  Scalar  bar,  3p.m.  (d-f)  Raman  and  PL  peak  position  mappir^  at  the  same  vibrating  modes  Scalar  bar, 
3tim.  Cg)  PL  spectra  recorded  from  uniform  strain  (black)  arrd  strain^vibratbn  locations  (red).  Inset  shows  a  "-8-meV  blue  shift  corresponding  to  a 
reduced  tensile  strain  (0.2%  less)  comparing  with  the  location  of  uniform  strain,  (h)  DFTcalculations  of  bandgap  shift  of  monolayer  MoS^  versus  tensile 
strain.  The  dots  represent  energy  shift  under  an  increasing  strain  from  DFT  calculations  and  the  trend  is  shown  by  the  dash  lir>e. 


PD  MS  (Fig,  3c,  under  4.B%  tensile  strain).  The  PL  intensity^ 
(Fig.  3d)  and  peak  position  (Fig.  3e)  mapping  were  obtained  for 
this  sample  under  an  increasing  tensile  strain  from  0  to  4.8%. 
Individual  PL  spectra  are  shown  in  Fig.  3£  The  corresponding 
energy  shifr  and  intensity  versns  tensile  strain  are  plotted  in 
Fig.  3g,  It  is  worth  noting  that  no  hysteresis  was  observed  for  the 
loading  and  unloading  pmcessesj  whidi  exclude  the  possibility  of 
sliding  effect  between  M0S2  and  substrates  (Supplementary 
Fig.  4). 

It  is  apparent  that  not  all  of  the  applied  strain  to  the  PDMS 
substrate  was  transferred  to  M0S2  crystals,  if  one  were  to  compare 
the  results  from  Figs  2g>h  and  3f,g.  The  heterogeneity  of  PL  at 
zero  strain  also  suggests  that  the  transfer  process  and  PDMS 
surface  properties  result  in  a  non  uniform  contact  between  the 
substrate  and  material.  To  better  explain  the  experimental  results 
and  estimate  the  strain  transferred  to  M0S2  layer  under  such 
tensile  tests,  a  3D  FE  model  was  developed,  as  sdiematicatty 
shown  in  Fig.  4,  Due  to  the  symmetry,  one  half  of  the  triangular 
sample  was  used  in  the  modd,  and  the  model  consists  of  one 
monolayer  M0S2  with  a  thickness  of  0.7  nm,  placed  on  the  surface 
of  PDMS  substrate  with  a  thickness  of  100  |im  (Fig.  4a  c). 

Figure  4d  h  shows  the  strain  distribution  in  the  triangular 
M0S2  different  applied  strains  on  PDMS.  The  tensile 


strain  is  dearly  transferred  from  the  substrate  to  the  M0S2  layer. 
The  strain  transferred  to  M0S2  layer  initiates  from  the  top  of  the 
triangle  (point  A  in  Fig,  4c)  and  then  propagates  down  along 
the  trian^e  and  finally  extends  to  the  bottom.  Interestingly,  the 
predicted  strain  distribution  across  the  M0S2  is  in  good 
agreement  with  the  PL  peak  position  mapping  shown  in 
Fig  3e,  indicating  that  the  PL  peak  position  distribution  is 
closely  rorrekted  to  the  local  strain.  The  strains  in  M0S2  layer  as  a 
function  of  the  applied  strain  on  PDMS  are  shown  in  Fig  4i, 
where  a  linear  relationship  is  predicted  within  the  ran^  of  5%  at 
which  the  substrates  have  still  not  yielded  under  such  applied 
strain  levels,  When  the  applied  strain  reaches  4.8%,  as  in  the 
experiment,  the  strain  in  M0S2  layer  varies  ranging  from  0,50, 
Q27  and  0. 14%  at  points  A,  B  and  C  (Fig  4c),  respectively.  Thus, 
up  to  10%  of  the  applied  strain  has  been  effectively  transferred  to 
M0S2  layer  in  die  case  of  assuming  perfect  van  der  Waals  bonding 
between  the  M0S2  and  the  substrate.  These  results  correlate  well 
with  the  eailier  PL  peak  position  mapping  of  the  CVD  grown 
M0S2  crystal  with  strain  variations  and  corresponding  DFT 
calculations,  confirming  the  validity  of  the  proposed  method  for 
strain  detection  and  bandgap  engineering  evaluation  using  the  PL 
peak  position  mapping  technique.  They  also  demonstrate  the  role 
of  substrate  and  crysti  geometry. 
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Figure  3  |  Controlled  tensile  strain  in  single-crystal  CVD  l\/loS2.  (a)  Schematic  of  PDMS/M0S2  stacks  for  tensile  strain  experiment.  The  uniaxial 
tensile  force  is  applied  along  the  direction  of  the  black  arrows,  as  indicated.  (b,c)  Optical  images  of  M0S2  without  and  under  4.8%  strain,  respectively. 
Black  arrows  show  the  direction  of  the  applied  force.  Scalar  bar,  6  pm.  (d,e)  PL  intensity  and  peak  position  mapping  under  an  increasing  tensile  strain  from 
0  to  4.8%.  Scalar  bar,  4  pm.  (f)  Representative  PL  spectra  under  each  strain,  collected  from  the  centre  of  the  triangle,  (g)  Energy  shift  (red)  and 
intensity  (blue)  of  the  PL  peak  (strongest  PL  peak  from  the  A  exciton)  versus  the  tensile  strain  of  PDMS.  Multiple  spectra  were  recorded  at  each  strain. 


In  addition,  it  is  important  to  calculate  the  effective  strain 
transfer  from  the  polymer  to  M0S2  and  demonstrate  how  to 
promote  the  strain  transfer  for  bandgap  engineering.  In  our 
experiments,  with  4.8%  tensile  strain  applied  to  the  PDMS 
substrate,  only  a  15  meV  shift  was  observed  in  the  PL  spectra 
near  the  centre  of  the  triangle,  corresponding  to  ~0.5%  tensile 
strain  according  to  DPT  calculations.  It  should  be  noted  here  that 
DFT  LDA  (local  density  approximation)  usually  produces  a 
relatively  lower  value  than  the  true  bandgap,  yet  the  change  in 
bandgap  can  be  more  accurately  determined,  although  the 
absolute  value  is  underestimated^^.  Our  work  suggests  that  only 
~  10%  of  the  applied  strain  is  effectively  transferred  to  M0S2 
when  using  PDMS  as  a  substrate.  Note  that  this  strain  transfer 
efficiency  was  also  confirmed  by  the  earlier  FE  simulation  (Fig.  4). 
Therefore,  it  will  be  of  practical  importance  to  further  investigate 
the  strain  loss  mechanisms  in  such  a  system.  We  calculated  the 
strain  in  the  M0S2  layer  for  various  substrates  with  Young’s 
modulus  ranging  from  0.35  KPa  to  1  GPa.  The  maximum  transfer 
strain  versus  Young’s  modulus  is  shown  in  Fig.  4j.  If  the  Young’s 


modulus  of  a  substrate  is  below  1  MPa,  the  transferred  strain  is 
about  0.5%  for  an  applied  strain  of  5%.  With  the  increasing 
substrate  modulus,  the  transferred  strain  quickly  increases  to 
4.7%  for  the  same  applied  strain,  indicating  that  almost  90%  of 
the  applied  strain  is  transferred  to  the  M0S2  layer.  As  a 
controllable  experiment,  we  transferred  the  M0S2  to  polyvinyl 
alcohol  with  Young’s  modulus  around  650  MPa,  around  three 
orders  of  magnitude  higher  than  PDMS  ( ~  430  KPa).  It  turns  out 
~  60%  strain  that  has  been  effectively  transferred  to  M0S2  based 
on  a  linear  extrapolation  from  Fig.  2h  (Supplementary  Fig.  5). 
This  result  supports  our  FE  simulation  predictions.  Thus, 
increasing  the  substrate  modulus  could  significantly  increase  the 
strain  in  M0S2  layer  if  a  high  quality  interface  is  preserved.  This 
result  provides  a  very  useful  guideline  for  the  fabrication  of  highly 
tunable  M0S2  electronic  devices. 


Structural  heterogeneity  in  CVD-grown  M0S2  monolayers. 

In  addition  to  strain  heterogeneity,  GBs  are  one  of  the  most 
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Figure  4  |  FE  simulation  of  tensile-strained  MoS2>  (a)  Schematic  of  the  tensile  test  sample,  (b)  FE  model  of  the  tensile  test  sample  consistirig  of 
a  lOO-jim-thick  substrate  (PDMS)  and  (c)  a  monolayer  M0S2  (0.7  nm)  on  the  PDMS  substrates.  Half  of  the  tensile  test  sample  was  modelled  due  to 
symmetry.  (d-h>  Snapshots  of  FEA  simulation  of  MoS^  tensile  tests:  applied  strain  (d)  0%,  (e)  15%,  (f)  2.2%,  (g)  3.4%  and  (h>  4.8%.  (i>  Strain  in 
M0S2  layer  as  a  function  of  the  applied  strain  on  the  substrate  at  points  A  (black),  B  (red)  and  C  (blue),  (j)  Maximum  strain  in  M0S2  layer  as  a  fuiKtion 
of  substrate's  Young's  Modulus  for  25%  (black)  and  5%  (red)  tensile  strain. 


significant  structural  features  in  CVD  grown  polycrystalline 
M0S2  samples  as  shown  in  Fig.  lb  d,  Attanpts  to  gro  w  large  area 
crystals  of  M0S2  has  rendered  an  unavoidable  polycrystallinity 
that  caBs  for  better  and  simi^er  identification  m^od  for 
detecting  the  GBs  and  understanding  their  role  in  the  electronic 
and  optical  properties  of  the  material.  The  experimental 
and  theoretical  work  on  the  GBs  of  M0S2  has  revealed  a  rich 
diversity  in  their  structure^^’^^.  It  is  known  that  monolayer 
M0S2  has  a  hexagonal  lattice  structure  similar  to  graphene 
and  h  BN.  However,  in  contrast  to  single  atom  thick  graphene, 
M0S2  nio  nolayers  consist  of  two  layers  of  sulfur  atoms 
sandwiching  a  layer  of  molybdenum  atoms.  This  results  in  a 
variety  of  dislocation  cores^  constituting  the  grain  boundary 
configurations  in  MoS2»  as  exposed  to  the  most  common  5  7  ring 
configuration  in  graphene^  In  addition,  variations  in  the  lattice 
orientation  mismatch  at  the  boundaries  and  non  bonding 
second  layer  growth  are  common  futures  observed  in  C]IVD 
grown  large  area  M0S2  films'^.  It  was  also  found  that  the  twin  GB 
in  M0S2  may  exhibit  metallic  behaviour^^.  The  most  reliable 
method  for  identification  of  GBs  has  traditionally  beoi  the 
hi^  resolution  TEM  techniques,  which  can  directly  inspect  its 
detailed  atomic  structural  characteristics.  This  approach  provides 
a  unique  atomic  level  perspective  of  the  M0S2  boundaries,  but  is 


not  suitable  for  direct  examination  of  large  area  samples/devices 
that  are  ready  for  subsequent  optical  and  electrical  measurements. 
(Consequently,  optical  based  microscopy  techniques  for 
identification  and  characterization  of  GBs  have  been  explored 
using  hi^  resolution  PL  mapping. 

Non  contact  optical  techniques  such  as  Raman  ^lectro&copy 
can  be  employed  to  map  the  GBs  of  graphene.  Since  the  Raman 
spectra  for  M0S2  is  v^  sensitive  to  the  number  of  layers, 
it  has  been  proposed  as  a  suitable  approach  for  detection  of 
over  layered  M0S2  GBs  (Supplementary  Rg.  6).  Nevertheless* 
for  monolayer  M0S2*  this  tedmique  is  not  efficient  in  detecting 
GBs  because  the  Raman  signal  is  very  weak  and  stron^y 
dependent  on  the  laser  power  for  monolayer  M0S2.  PL  variations 
at  the  GB  have  also  been  exploredi  however,  little  is  known 
about  the  fundamental  reasons  for  the  observed  changes. 
As  discussed  earlier,  the  PL  peak  position  mapping  is  proven  to 
be  an  effective  method  to  detect  intrinsic  strains  in  monolayer 
M0S2  (ref.  21),  It  is  wdl  known  that  the  grain  boundary  in  2D 
material  consists  of  lined  up  dislocations,  which  naturally  may 
introduce  localized  strain  in  their  vicinity.  In  addition,  the  variety 
in  CVD  grown  crystal  geometries  and  grain  boundary  types  can 
be  utilized  to  further  gather  insight  into  their  role  in  M0S2  PL 
characteristics. 
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The  edges  and  twin  GBs  have  a  higher  intensity  (light  yellow),  while  rest  are  depicted  in  orange,  (c)  PL  peak  position  mapping.  PL  spectra  are  collected 
at  the  location  marked  five  dots  with  different  colours,  (d)  PL  spectra  collected  from  points  indicated  in  c.  Each  spectrum  corresponds  to  the  location 
in  c  with  the  same  color.  (e,f)  PL  intensity  and  peak  position  line  profiles  as  indicated  by  the  white  dash  lines  crossing  two  GBs  and  two  edges. 


To  analyse  the  PL  response  of  GBs  in  polycrystalline  M0S2,  we 
examine  a  representative  six  pointed  M0S2  star  with  six  twin  GBs 
(Fig.  5a).  The  PL  intensity  mapping  (Fig.  5b)  reveals  the  enhanced 
intensity  characteristics  of  six  mirror  twin  GBs.  Such  twin 
structures  have  also  been  confirmed  by  the  PL  peak  position 
mapping  (Fig.  5c).  The  same  structure  has  been  observed  under 
dark  field  TEM^.  It  was  found  that  the  outer  edges  are  oriented 
along  the  Mo  z  z  (molybdenum  zigzag)  direction,  while  GBs 
occur  along  the  S  z  z  (sulfur  zigzag)  direction  (Supplementary 
Fig.  7)^.  Nevertheless,  these  GBs  represent  unique  defect 
characteristics.  PL  line  profiles  and  five  individual  PL  spectrums 
along  the  lines  indicated  in  their  respective  scans  are  recorded 
(Fig.  5d  f).  Their  locations  are  labelled  in  Fig.  5c.  Three  of 
the  points  are  located  at  the  pristine  M0S2  and  two  of  them  at  the 
GBs.  An  enhancement  in  the  PL  up  to  250%  is  observed  at  the 
GBs  (Fig.  5d,e).  The  highest  PL  intensity  is  found  at  the  GBs  and 
edges,  which  is  opposite  to  the  report  from  Van  der  Zande  et  al? 
who  observed  a  lower  PL  intensity  along  the  GBs  in  monolayer 
M0S2.  The  physical  scenario  of  PL  along  the  GBs  is  very 
complicated.  According  to  previous  studies,  defects  at  the  edges 
and  the  GBs  may  contribute  the  enhancement  of  PL  intensity. 
New  sub  bandgap  emission  peak  as  well  as  increase  in  overall  PL 
intensity  were  found  as  a  result  of  the  vacancy  generation^ ^  via  a 
defect  assisted  process^^.  Here,  our  results  also  imply  that  a  local 
interaction  also  plays  an  important  role  in  the  PL  variations.  We 
propose  that  the  non  uniform  strain  due  to  the  thermal  mismatch 
of  the  substrate,  randomness  of  interface  properties  and  strain 
transfer  and  tensile  straining  native  to  GBs  combined  make  the 
most  significant  contribution  to  the  optical  bandgap  shifts.  We 
find  that,  in  Fig.  5,  the  PL  peak  shift  is  ~  5  meV  at  the  GBs  and 
~  15meV  at  the  edges  in  monolayer  M0S2.  This  small  blue  shift 
indicates  that  no  matter  what  type  of  defects  exist  at  the  GBs, 
their  behaviour  tends  to  be  very  similar.  This  further  supports  our 


hypothesis  that  strain  at  the  GBs  is  possibly  the  main  reason  for 
the  observation  of  changes  in  the  PL  characteristics  of  M0S2.  To 
confirm  the  strain  in  GBs,  we  collected  low  temperature  PL 
spectra  at  non  GB  and  GB  locations.  The  shifts  have  been  found 
for  both  locations  and  the  difference  in  peak  location  increases, 
demonstrating  that  the  PL  spectra  at  the  GBs  are  more  sensitive 
to  the  change  of  strain  (Supplementary  Fig.  8).  We  have 
examined  the  temperature  dependency  of  the  bandgap  and  PL 
positions  of  M0S2  atomic  layers  using  a  modified  Varshni  model 
(Supplementary  Fig.  8)^^.  These  results  show  a  good  fit  for  the 
measurement  of  the  optical  gap  in  the  crystalline  region  and  the 
grain  boundary  regions  of  M0S2,  implying  that  the  changes  seen 
in  the  PL  characteristics  are  not  results  of  disorder  and 
localizations  in  the  excitonic  transitions.  Also,  not  only  for  the 
symmetric  M0S2  stars,  similar  results  have  been  seen  from  the 
asymmetric  polycrystalline  M0S2  flakes  with  GBs  (Supplementary 
Fig.  9).  In  addition,  we  have  also  noticed  that,  in  some  M0S2 
polycrystalline  flakes,  the  GBs  were  broken  induced  by  tensile 
strain.  It  further  confirms  the  vulnerability  of  GBs  in  M0S2 
(Supplementary  Fig.  10)  and  the  impact  of  strain. 

Discussion 

In  this  work,  we  have  used  PL  mapping  to  study  strain  and 
structural  heterogeneity  in  PL  characteristics  of  CVD  grown 
M0S2  monolayers.  We  observe  and  quantify  the  native  non 
uniform  strain  in  such  samples  and  propose  that  the  difference  in 
thermal  expansion  coefficient  of  M0S2  and  the  supporting  Si02 
substrate  leads  to  an  intrinsic  tensile  strain  in  M0S2  during  the 
fast  cooling  process  at  the  final  stage  of  the  growth.  A  relationship 
between  the  PL  peak  position  shifts  and  strain  is  established  and 
can  be  used  for  strain  detection  and  bandgap  engineering 
evaluation.  Further  examination  of  substrate  mediated  straining 
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process  provides  important  insight  into  the  role  of  substrate  in 
effective  strain  transfer.  Through  FE  simulation  and  controlled 
tensile  experiments,  we  find  that  for  PDMS  substrates  only  ~  10% 
of  the  applied  strain  can  be  effectively  transferred  to  M0S2. 
Moreover,  we  demonstrate  how  strain  concentration  propagates 
in  a  triangular  M0S2  monolayer  crystal  and  propose  substrates’ 
Young’s  modulus  as  dominating  mechanisms  of  strain  loss.  A 
substrate  with  a  modulus  >  1  GPa  may  effectively  transfer  almost 
90%  of  the  applied  strain  to  monolayer  M0S2  when  perfect 
interfacial  contact  conditions  are  provided.  In  addition,  we 
examine  the  discrepancies  in  the  literature  about  the  role  of  GBs 
in  PL  characteristics  of  M0S2.  Our  results  suggest  that  thermally 
induced  non  uniform  tensile  strain  plays  a  significant  role  in 
grain  boundary  optical  properties.  Our  results  also  provide  a 
promising  way  for  fabrication  and  study  of  highly  tunable  M0S2 
electronic  devices  and  mechanisms  of  strain  transfer  in  2D  atomic 
layers. 

Methods 

CVD  of  monolayer  M0S2.  To  synthesize  M0S2,  the  precursors,  M0O3  films,  were 
prepared  through  the  filtration  of  their  highly  crystalline  nanoribbons.  These  rib¬ 
bons  were  produced  hydrothermally  through  a  process  widely  used  for  synthesis  of 
this  type  of  single- crystalline  material.  Sodium  molybdate  (Na2Mo04)  or 
hexaammonium  heptamolybdate  tetrahydrate  (NH4)6Mo7024  •  4H2O  was  used  as 
precursors  for  the  preparation  of  M0O3  nanoribbons.  We  dissolve  1.2  g  of  these 
precursors  in  nitric  acid  and  transfer  to  a  Teflon  container  and  heat  it  at  170  °C  for 
1  2  h.  These  M0O3  nanoribbons  have  a  high  aspect  ratio,  roughly  20  pm  in  length, 
1  2  pm  in  width  and  thickness  in  the  range  of  10  40  nm.  After  the  preparation  of 
these  M0O3  precursors,  the  ribbons  were  Altered,  and  large-area  Aims  were  formed, 
cut  into  pieces  and  dispersed  on  silicon  substrates.  Several  bare  and  clean  sub¬ 
strates,  designate  for  the  growth  of  M0S2,  and  the  M0O3  Aim  were  placed  close  to 
each  other  at  the  centre  of  the  furnace  vented  with  nitrogen  at  200  seem.  A  con¬ 
tainer  was  Ailed  with  0.8  1.2  g  of  sublimated  sulfur  and  placed  at  the  opening  of  the 
furnace  at  a  location  that  reaches  an  approximately  maximum  temperature  of 
600  °C.  The  heating  procedure  corresponding  to  the  temperature  at  the  centre  of 
the  furnace  starts  with  a  gradual  heating  from  room  temperature  to  550  °C  in 
30  min  at  a  rate  of  ~20°Cmin“  ^  Then,  at  a  slower  pace,  ~5°Cmin“  ^  the 
chamber  is  heated  to  850  °C.  AAer  10  15  min  at  this  temperature,  the  furnace  is 
cooled  naturally  back  to  room  temperature.  At  temperatures  close  to  550  °C  at  the 
centre,  sulfur  will  slowly  start  to  evaporate  and  the  central  interaction  between  this 
element  and  M0O3  follows. 

Both  Raman  and  PL  data  were  collected  using  a  532-nm  laser  excitation  focused 
through  a  x  100  objective  lens.  This  resulted  in  a  spot  size  of  ~342nm  where 
scanning  was  performed  using  200  nm  steps  on  a  piezoelectric  stage  to  obtain  high- 
resolution  images.  All  Raman  and  PL  spectra  presented  in  this  manuscript  for 
M0S2  on  Si02/Si  were  taken  at  an  incident  laser  power  of  12.5  pW,  which  was 
sufflciently  low  to  avoid  any  shiAing  of  the  Raman  modes  during  scanning. 
Mapping  on  the  PDMS  substrate  was  performed  at  an  incident  laser  power  of 
50  pW  to  improve  the  signal  to  noise  ratio.  The  PL  peak  position  was  found  to  be 
much  less  sensitive  to  the  laser  power  density.  SpeciAcally,  the  laser  power 
dependence  on  Raman  peak  position  and  intensity  is  shown  in  Supplementary 
Figs  11  and  12  and  shows  that  the  layer  number  can  easily  be  incorrectly 
determined  when  scanning  at  high  power  densities. 


DFT  calculations.  DFT  calculations  within  the  LDA  were  performed  using 
the  projector-augmented  wave  (PAW)  method^^’^^,  as  implemented  in  Vienna 
Ab-initio  Simulation  Package^^’^F  Twelve  electrons  were  included  for  Mo 
pseudopotential.  A  vacuum  layer  >  12  A  was  introduced  to  make  the  spurious 
interactions  negligible.  Allxllxl  Monkhorst-Pack  k-point  mesh  centred  at  the 
r -point  was  adopted  for  the  Brillouin  zone  integration.  With  an  energy  cutoff  of 
280  eV,  a  plane-wave-based  total  energy  minimization  scheme^^  was  employed 
until  the  force  on  each  atom  was  <0.01  eV  A“  ^  The  uniaxial  strain  was  applied 
either  in  the  armchair  (AC)  or  zigzag  (ZZ)  direction. 


3D  FEA.  The  discretized  FE  models  contain  about  53,254  nodes  and  49,584 
elements  with  very  Ane  meshes  in  the  contact  region.  The  Anal  mesh  density 
was  determined  through  a  series  of  convergence  studies.  Appropriate  boundary 
conditions  were  used  along  two  edges  of  the  specimen  to  simulate  the  loading 
strain  during  testing.  The  interface  between  M0S2  and  PDMS  were  modelled  using 
perfect  bonding.  The  calculations  were  performed  using  the  commercial  FE 
package  ABAQUS  (version  6.10).  Young’s  modulus  and  Poisson’s  ratio  for 
PDMS  are  350  KPa  and  0.5,  and  300  GPa  and  0.3,  respectively,  were  used  in 
the  simulation. 
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